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ABSTRACT 

High precision fast photometry from ground-based observatories is a challenge due 
to intensity fluctuations (scintillation) produced by the Earth's atmosphere. Here we 
describe a method to reduce the effects of scintillation by a combination of pupil 
reconjugation and calibration using a comparison star. Because scintillation is pro- 
duced by high altitude turbulence, the range of angles over which the scintillation is 
correlated is small and therefore simple correction by a comparison star is normally 
impossible. We propose reconjugating the telescope pupil to a high dominant layer 
of turbulence, then apodizing it before calibration with a comparison star. We find 
by simulation that given a simple atmosphere with a single high altitude turbulent 
layer and a strong surface layer a reduction in the intensity variance by a factor of 
~ 30 is possible. Given a more realistic atmosphere as measured by SCIDAR at San 
Pedro Martir we find that on a night with a strong high altitude layer we can expect 
the median variance to be reduced by a factor of ~ 11. By reducing the scintillation 
noise we will be able to detect much smaller changes in brightness. If we assume a 
2 m telescope and an exposure time of 30 seconds a reduction in the scintillation noise 
from 0.78 mmag to 0.21 mmag is possible, which will enable the routine detection of, 
for example, the secondary transits of extrasolar planets from the ground. 

Key words: atmospheric effects - techniques: photometric 



1 INTRODUCTION 

High precision fast photometry is key to several branches of 
research including (but n ot limited to) the st udy of extra- 
solar planet transits (e.g. Charbonneau et al.l 2000), stellar 
seismology ([Christensen-Dalsgaard et al.ll2006l) and the de- 
tection of small Kuiper belt objects (e.g. ISchlichting et al.l 
2009). The difficulty with such observations is that, although 
the targets are often bright, the variation one wishes to de- 
tect is often very small (typically millimagnitudes or less) 
and hence the signal to noise ratio is not limited by the 
detector or sky but by intensity fluctuations (scintillation) 
produced by the Earth's atmosphere. For this reason fast 
photometers are generally put in space (e.g. CoRoT, Kepler 
and PLATO). 

Extrasolar planetary transits can be detected from the 
ground. However the measurement of the secondary eclipse 
(i.e. where the planet goes behind the star) is a challenge. 
Such observations are crucial, as only the secondary eclipse 
can give information on the planetary atmosphere, includ- 



ing the temperature and albedo ([Knutson et al.|[2007h . Sec- 
ondary eclipses were detected fo r the first time from spac e 
in 2005 using Spitzer at 3 /xm (|Charbonneau et alJliooEl ). 
There has been a great deal of effort to detect secondary 
eclipses from the ground, but for years no detections were 
made (in large part due to scintillation noise). Finally, in 
2009, the first ground-based detections were made, but these 
relied on near-IR measurements and had to target the most 
bloated, closest (t o their host star) exoplanets to maximise 
the eclipse signal (Sin g &; Lop ez- Morales! 120091 ) . Since then 
a few other exoplanets have had secondar y eclipses detected 
from t he ground in this way. As noted bv lDeming &; Seagerl 
(2009), secondary eclipses recorded in visible light in addi- 
tion to IR measurements are crucial if we are to understand 
the relative contribution of thermal emission and reflected 
light, and the planetary albedo. 

Time averaging the intensity will reduce the scintillation 
noise by an a mount proportional to the square root of the 
exposure time ([Dravins et al.|[T998l ) , but this will often result 
in saturating the CCD which then requires de-focusing the 
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telescope to distribute the image of the star over more pixels. 
De- focusing has certain advantages, such as reducing the im- 
pact of pixel-to-pixel and intra-pixel sensitivity variations, 
but it also significantly increas es the sky and CCD readout 
noise (|Southworth et all l2009h . In addition, de- focusing is 
not routinely possible on some telescopes (e.g the VLT) and 
it can not be done with crowded fields. More importantly 
for fast photometry, time averaging can also only be used 
in circumstances where the intrinsic variability of the target 
has a much longer time scale than the scintillation. As scin- 
tillation is caused by the spatial intensity fluctuations cross- 
ing the pupil boundary, the time sc ale is determined by t he 
wind speed of the turbulent layer. iDravins et al ] (|l997al lbl 
1998) studied the temporal autocorrelation of the scintilla- 
tion pattern at astronomical sites and found that the power 
is mainly located between 10 and 100 Hz but actually spans 
many orders of magnitude. 

Differential photometric measurements ca n be made by 
normalising with a nearby comparison star fe.g. lHenrv et all 
2000). This is not to reduce the scintillation but to correct 
for transparency variations in the atmosphere. However, this 
actually makes the scintillation noise worse as it is inherently 
caused by high altitude layers and therefore will have a very 
small angle of coherence in the optical (typically ~ l"). Here 
we propose a technique, called "conjugate-plane photome- 
try" to reduce scintillation noise by increasing the angle of 
coherence up to ~ 0.5°, allowing the intensity variations of 
the target star to be corrected by a comparison star. Our 
technique offers a relatively simple way of routinely obtain- 
ing space-quality photometry from the ground for a fraction 
of the price and with much larger telescope apertures. 

In section 2 we describe the scintillation reduction 
method. Section 3 shows the results of simulations of our 
correction technique. The expected performance of the sys- 
tem for a theoretical extrasolar planet transit, and simula- 
tion results using a real atmospheric profile are shown in 
section 4. Finally in section 5 we discuss the design of a 
prototype which will be tested at the NOT on La Palma in 
September 2010. 



2 SCINTILLATION CALIBRATION 

High altitude turbulence in the atmosphere distorts the 
plane wavefronts of light from a star which is effectively 
at infinity. As the wavefronts propagate these phase aber- 
rations evolve into intensity variations which we view with 
the naked eye as twinkling. Wavefronts incident on a tele- 
scope pupil have both phase variations, caused by the in- 
tegrated effect of light passing though the whole vertical 
depth of the atmosphere, and intensity variations, caused 
predominantly by the light diffracting through high altitude 
turbulence and interfering at the ground. Phase variations 
are normally considered more significant as they dramati- 
cally affect the spatial resolution of images, and this has led 
to the development of adaptive optics. The intensity varia- 
tions across the pupil are effectively averaged together when 
the light is focused and therefore have less effect. A larger 
aperture implies more spatial averaging (which is why stars 
twinkle less when observed through a telescope than with 
the naked eye). However, these small intensity fluctuations 
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Figure 1. A spherical wavefront from a star will appear flat as 
it enters the the atmosphere. In the absence of turbulence this 
flat wavefront will be collected by the telescope pupil (left). In 
the presence of turbulence the wavefront will diffract through 
the refractive index variations which accompany the turbulent 
motion in the atmosphere. The wavefront will then interfere with 
itself at the ground and cause intensity fluctuations. A simplified 
geometrical model is shown on the right. The scintillation noise 
occurs when extra light is focused into the telescope pupil or when 
light is focused away from the pupil by the turbulent atmosphere. 

do become significant when one is concerned with high pre- 
cision photometry. 

Consider now the effect of these intensity variations in 
more detail. If we ignore diffraction, then a flat wavefront 
which is the same size as the telescope pupil at a given 
high altitude, in the absence of atmospheric turbulence, will 
propagate in a direction normal to the wavefront and will 
be collected by the telescope pupil (see figure [1]) . Now con- 
sider the effect of atmospheric distortion. Phase aberrations 
cause diffraction in different directions and hence produce 
scintillation. Effectively light from one part of the original 
wavefront is redirected to other parts of the pupil. This in 
itself is not a significant problem for photometry, as the in- 
tegrated intensity across the pupil is the same. The problem 
occurs either when rays from the wavefront at high altitude 
propagate away from the telescope pupil, and are lost, or 
conversely when high altitude areas away from the telescope 
pupil area propagate into the telescope pupil at the ground. 
These effects lead to a decrease and increase in intensity, re- 
spectively, and at any one instant both of these effects will 
be occurring (see figure [1]) . The turbulence is blown across 
the field of view of the telescope producing an overall change 
in intensity as a function of time. 

As a thought experiment, to show the basic concept be- 
hind our proposal, if we could place an aperture which is 
smaller than the telescope pupil in the sky at the altitude 
of high turbulence then this change in intensity could be 
dramatically reduced. In this case, the rays that would have 
been deflected away from the area of the pupil would still 
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Figure 2. By placing an aperture at the altitude of the turbu- 
lent layer we can reduce the scintillation noise. It will now be 
impossible for any light from outside of the telescope pupil to be 
focused into the collecting area. It will also be unlikely for any 
parts of the wavefront to be focused off-axis by such a degree as 
to escape from the collecting area all together. These two situa- 
tions are shown in red. These rays - which would normally be the 
ones producing a change in the overall intergrated intensity - are 
blocked. 



be collected by the (larger) telescope pupil, and as the an- 
gle of diffraction is small no rays would be deflected into 
the telescope pupil because of the aperture (see figure [2]). 
Increasing the size difference between the aperture in the 
sky and the telescope pupil would improve the scintillation 
rejection, but would also lead to increased loss of signal, and 
clearly a balance between the two effects would need to be 
found. 

Placing an aperture at a high altitude in the sky is 
clearly an impractical proposal, but we can produce a sim- 
ilar effect using optics after the telescope focus. Figure [3] 
shows how reconjugation can be produced by observing the 
beam in a different plane downstream from the telescope fo- 
cus. The high altitude turbulent layer is reimaged onto an 
aperture which is slightly smaller than the equivalent size of 
the full telescope pupil. Consider again the simplified case 
of a single layer of turbulence at a high altitude. As already 
described, this produces scintillation in the entrance pupil 
of the telescope. If we reimage the high altitude layer at a 
conjugate plane then the rays will have propagated so as 
to "undo" the scintillatio n and we would v iew an approxi- 
mately uniform intensity (|Fuchs et al.lll998h . High altitude 
areas of the wavefront, which in the absence of turbulence 
would fall outside of the telescope pupil, can be diffracted by 
the turbulence and interfere to cause intense regions within 
the pupil area. This light would image in the conjugate 
plane outside of the aperture and can be easily rejected by 
the mask. High altitude areas of the wavefront which are 
diffracted by the turbulence and interfere to cause intense 
areas at the ground outside of the telescope pupil are lost 
and will show up as areas of decreased intensity towards the 
edge of the reimaged wavefront. This effect can also be re- 
jected with a mask at the reimaged altitude which is slightly 
smaller than the pupil size. The remaining light within this 
mask will be approximately of uniform intensity and scintil- 
lation free. 

The above description has ignored two important ef- 
fects, namely diffraction and turbulence from other atmo- 
spheric layers (predominantly low altitude turbulence). As 
well as high altitude turbule nce most astronomical sites will 
have a strong surface layer (|Osborn et alJl2Ql(il ; IChun et al.l 
2009) and possibly turbulence at intermediate altitudes as 
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Figure 3. Ray diagrams for conjugation positions. The black 
lines show the rays for an object at infinity. The top diagram 
shows the conjugate position of the telescope pupil. Every point 
in this plane will be an image of a point on the telescope pupil 
(as shown by the red lines). The lower diagram shows that by 
moving the observation plane towards the collimating lens then 
an image of the wavefront at a height h above the telescope will 
be produced. If a camera is in a position such that it is in the 
image plane of the turbulent layer it is at the conjugate altitude 
of that layer. In practice subsidiary optics may also be used, but 
this diagram shows the principle. 
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Figure 4. In differential photometry the intensity of the target 
star is calibrated by the intensity from a second comparison star. 
As the scintillation is caused by high altitude turbulence the two 
light cylinders do not sample the same turbulence and hence there 
will be very little correlation between the two. By conjugating the 
telescope to the high-altitude layer we remove the scintillation 
from this layer and it is replaced by scintillation from the surface 
turbulent layer instead. However, as the two light cylinders sample 
the same region of turbulence near the ground they will have very 
similar scintillation patterns, allowing one to be corrected by the 
other. The angle of separation of the two stars can be large as the 
surface layer is generally found to be thin. 
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well. If we conjugate our system to the altitude of a high 
turbulent layer we will still see scintillation from other lay- 
ers. We will have effectively swapped scintillation caused by 
high altitude turbulence with scintillat i on ca used by turbu- 
lence close to the ground. iFuchs et al.l (| 19981 ) demonstrated 
that if a turbulent layer is below the conjugate plane (the 
surface layer for example) then a virtual reverse propaga- 
tion occurs over a distance z — \h — zq\, where zo is the 
conjugate altitude and h is the altitude of the turbulent 
layer. Therefore the surface layer will now cause scintilla- 
tion in the conjugate plane as it will have effectively prop- 
agated a distance zo. However a comparison star can be 
used to reduce the scintillation from the surface layer as 
they will both sample the same turbulent area, as shown in 
figure 21 This layer must also be quite thin to ensure the 
wavefront samples the same turbulence, and studies have 
demonstrated t hat this is th e case (it is often only a few 
IP's of met ers, lOsborn et all 2010, iTokovinin et all 2010, 
IChun et all 2009) meaning that the coherence angle is now 
very large (up to 0.5°). 

Figure [5] shows the effect of reconjugation of a sin- 
gle high altitude layer, including the effects of diffraction 
caused by the telescope pupil. The simulation assumed a 
single high altitude turbulent layer at 10 km with f C^dh = 
353 x 10 -15 m 1 / 3 , where C% is the refractive index structure 
constant and f C^dh is the integrated turbulence strength 
of the atmospheric layer. This corresponds to ro = 0.15 m, 
where ro is the Fried parameter and is a measure of the in- 
tegrated strength of the turbulence. It can be seen that the 
variations in intensity due to scintillation largely disappear 
in the reconjugated image of the high altitude layer - but 
that diffraction can clearly be seen. The diffraction rings 
are not completely circular as a result of the phase distor- 
tions in the wavefront at the telescope pupil. Figure \6\ shows 
simulated images of the reconjugated pupils at 10 km for a 
two-layer atmosphere (0 and 10 km) for two stars separated 
by 40" . The two images are very similar indicating that one 
may be used to calibrate the other. They are not identical, 
however, as they are being not being illuminated by a flat, 
uniform wavefront due to the high altitude turbulence (and 
not the finite thickness of the layer), and this introduces a 
source of error - as highlighted in the next section. 



3 THEORY AND SIMULATION RESULTS 

Assuming a single turbulent layer at 10 km and no other 
turbulence the wavefunction, ^, at the telescope pupil is 
given by, 

(x, y) = [K(z = +10 km) <g> exp (ifao)] P(x, y), (1) 

where z is the propagation distance, x and y are spatial 
co-ordinates, P(x,y) is the telescope pupil function, (j>h is 
the turbulent phase screen at altitude h km, denotes a 
convolution and K is the Fresnel propagation kernel, given 
by, 

K= j- exp (ikz) exp ^ [(x-x') 2 + (y - y') 2 ^J , (2) 

where k is the wavenumber A is the wavelength of the light 
and x and y and spatial co-ordinates in the observation 
plane located at a distance z. Positive z indicates a diverging 
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Figure 5. Simulated pupil intensity patterns at the telescope 
pupil (left) and at the conjugate altitude of the turbulent layer 
(right). The telescope pupil is 2.0 m in diameter and the turbulent 
layer has f C^dh = 353 x 10 -15 m 1 / 3 and is located at an altitude 
of 10 km. The intensity pattern at the conjugate altitude shows 
that the spatial intensity fluctuations have been removed but have 
been replaced by diffraction rings concentrated around the edges 
that also permeate throughout the pupil. 
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Figure 6. Pupil images conjugate to 10 km for two stars sepa- 
rated by 40 // . The spatial intensity fluctuations are a combination 
of the scintillation pattern from the surface turbulent layer and 
the diffraction pattern of the telescope pupil (figure [5] right). The 
two images have very similar intensity patterns as they are both 
formed by the propagation of the same area of surface layer. 

spherical wavefront and negative z is a converging spherical 
wavefront or a negative propagation. Therefore, the wave- 
function in the conjugate plane, ^'(x',y'), is found by a 
further propagation of the wavefront by a negative distance, 

V'(x',y') = 

K(z = -10 km)®[[K(z = +10 km) <g> exp (ifao)] P(x, y)] . 

(3) 

In the case of an infinitely large pupil, ty'(x' \y') = *&(x,y) 
and the pupil amplitude is flat. Therefore, by placing the 
aperture at the conjugate altitude of the turbulent layer we 
can reduce the scintillation caused by that layer. However, 
with a real aperture the intensity profile at the conjugate 
plane is not flat because the wavefront diffracts through the 
telescope pupil and causes diffraction rings at the edge of 
the pupil image which are a function of the turbulent phase 
screen. If we include a ground layer, 0n, the Fresnel propa- 
gation equation becomes, 

y'{x',y) = K(z = -10 km) 

<g> [[K(z = +10 km) ® exp (z^io)] exp (i<j> )P(x, y)] . (4) 
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The surface layer and telescope pupil are multiplied into the 
wavefront before the final convolution. This is why these ef- 
fects can not be de-coupled from the higher turbulent layers 
and the wavefront in the conjugate plane will therefore de- 
pend on the high altitude phase aberrations as well as the 
surface layer and will be different for the target and com- 
parison stars. In addition to the diffraction these residual 
intensity variations will limit the effectiveness of the tech- 
nique. 

Our conjugate-plane photometry concept has been sim- 
ulated using a Fresnel propagation wave optics simulation 
using the theory stated above and randomly generated phase 
screens. Scintillation is often quantified by the scintillation 
index, <7g cint , which is defined as the normalised variance of 
intensity fluctuations, crg cint — ((/ — (7)) 2 )/(/) 2 , where / 
is the intensity o f the image and (I) denotes the time av- 
eraged intensity (|Dravins et al.l ll997aT ) . Figure [7] shows the 
scintillation index as a function of aperture size for a few 
example cases. The first case shows the theoretical maxi- 
mum reduction found by suspending the aperture in the sky 
above the telescope (solid line). This is entirely unfeasible 
but places a maximum limit on the reduction of the variance. 
The black dot-dashed line shows the scintillation variance 
for differential photometry with the aperture in the conju- 
gate plane. Diffraction through the pupil means that light 
is redistributed in the pupil. Therefore, simply blocking the 
outer regions of the pupil will no longer remove most of the 
extra light and will result in a higher scintillation variance. 
The small shoulder in the curve at approximately 0.07 m 
coincides with the radius of the first diffraction ring. The 
red dashed lines show the scintillation variance with a high 
altitude layer and a surface layer which varies in strength. In 
this case a comparison star is required to normalise the scin- 
tillation. The strength of the surface layer is selected so that 
the ratio of C 2 (10 km)d/i/C 2 (0 m)dh is equal to 1, 2 and 4. If 
the surface layer is weaker than the high turbulent layer the 
residual intensity fluctuations will be lower and so the resid- 
ual scintillation will also be lower. The maximum median 
variance reduction factor for C 2 (10 km)dh/Cn(0 km)dh = 1 
(i.e. equal strength), 2 and 4 is 17, 23 and 47, respectively 
and is found at i^aperture — D te \ ~ 0.1 m, for a simulated 
telescope diameter of 2 m. 

The amplitude of the first diffraction ring is substan- 
tially larger than any others (as seen in figure [5}. The opti- 
mum aperture size is therefore one which blocks this ring but 
none of the others. This will minimise the residual diffrac- 
tion and retain a large pupil area. The radius of the first 
diffraction ring in the very near field is given by the radius 
of the first Fresnel zone, tf = V~\z, in this case 0.07 m and 
is independent of telescope size. 

The reduction in scintillation noise can be clearly seen 
in figure [8] which shows the normalised light curve for a se- 
quence of 200 frames from a simulation assuming a constant 
source intensity. The black line shows the original light curve 
with a variance of 1.6 x 10 -4 due to scintillation. The red line 
is the light curve after scintillation reduction and has a vari- 
ance of 6.4 x 10 -6 , a reduction factor of 25. The variance is in 
units of normalised intensity, AI/I. The simulation assumes 
an atmosphere with two turbulent layers, one at the ground 
and one at 10 km, both with / Cldh = 353 x 10~ 15 m 1/3 , 
the telescope diameter was 2 m and there was no temporal 
averaging. 
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Figure 7. The solid line shows the scintillation variance as a 
function of aperture size for an aperture suspended in the sky 
10 km above a 2 m telescope. In this case it is possible to reduce 
the scintillation variance to effectively zero. The black dot-dashed 
line shows the scintillation variance for a single high-altitude tur- 
bulent layer with the aperture in the conjugate plane. The per- 
formance is not as good as the solid line due to the diffraction 
from the telescope pupil. The red dashed lines show the scintilla- 
tion variance for the aperture in the conjugate plane of the high 
turbulent layer and with a surface layer with strengths equal to 
C^(zo)dh, 2 x C^(zo)dh and 4 x C^(zo)dh : where zo is the conju- 
gate altitude, with C^(zo)dh = 3.5 x 10 _13 m 1 / 3 . The data points 
and error bars are the mean and standard errors of 20 simulations, 
each with unique and randomly generated phase screens. 

A mis-conjugation of the aperture will result in less than 
optimal performance. Figure [9] shows the factor by which the 
scintillation variance is reduced as a function of conjugate 
altitude for turbulent layers at m and 10 km. In this case 
the curve has a full width at half maximum of approximately 
3.5 km. This will be narrower for turbulent layers at lower 
altitudes and wider for higher altitudes. Knowledge of the 
contemporaneous turbulence profile is therefore essential to 
ensure that the aperture is conjugate to the correct altitude. 



4 PERFORMANCE ESTIMATION 

The Monte- Carlo simulations are useful to examine the per- 
formance for a particular parameter set. However, it is very 
inefficient for modelling the performance as a function of 
time for real turbulence profiles with many turbulent layers. 
To do this an analytical estimate of the intensity variance 
for a given parameter set is required. 

If the pupil is much larger than the Fresnel radius 
(D ^> VAzn) the intensity variance due to scintillation, 
ofr int, can be predicted us ing the theoretical model described 
bv lDravins et al.l (|l997bh . 

_ i r°° 

Cscint oc D tel 3 (sec7) 3 / Cl{h)h 2 dh, (5) 
J o 

where 7 is the zenith angle. The scintillation index is then 
independent of wavelength and proportional to the altitude 
of the turbulent layer squared and the strength of the tur- 
bulent layer. We can calculate the scintillation index due to 
all of the turbulent layers assuming the pupil is conjugate to 
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Figure 8. An example simulated light curve. The black line is 
the intensity pattern from a simulation observing a star with a 
2 m class telescope through the atmosphere with a turbulent layer 
at 10 km and m, both with / C^dh = 353 x 1CT 15 m 1 / 3 . The 
exposure time of each frame is short so that there is no temporal 
averaging of scintillation. The red line shows the scintillation cor- 
rected light curve. In this case the intensity variance is reduced 
from 1.6 x 1CT 4 to 6.4 x 1CT 6 , a factor of 25. The residual noise 
is due to the uncorrected scintillation. 
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Figure 10. SCIDAR turbulence profile, i.e. height above sea level 
against time, where the colour indicates the strength of the turbu- 
lence, from 2000 May 19 at San Pedro Martir. The profile shows 
a dominant layer at approximately 10 km throughout the night. 
San Pedro Martir is located at 2800 m above sea level. 
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Figure 9. Ratio of intensity variance for normal differential pho- 
tometry (cr^ iff ) and scintillation corrected photometry (c^ orr ) ver- 
sus conjugate altitude for an atmosphere with layers at m and 
10 km, both with / C^dh = 353 x 10 -15 m 1 / 3 and a telescope 
diameter of 2 m. The curve is Lorentzian with a FWHM of ap- 
proximately 3.5 km. At conjugate altitude m we measure an 
improvement in the intensity variance of ~0.5, i.e. the variance 
is actually increased. This is because the pupil size is reduced by 
the apodizing mask. The data points and error bars are the mean 
and standard deviation of 20 simulations, each with unique and 
randomly generated phase screens. 

an altitude, zq. In this case the scintillation index, a% Q1 at a 
given altitude can be calculated using a modification to the 
scintillation index equation (equation [5]) , 

Czq 

DJ (sec 7 )" / Ct (h) (h - zoY dh, (6) 

JSL 

where (h — zq) is the separation between the layer altitude 



and the conjugate altitude, ignoring the surface layer as this 
will be dealt with separately. 

The corrected residual scintillation variance, cr 2 ovv , will 
be dominated by this but we also add noise terms due to the 
pupil diffraction and the surface layer. These noise sources 
are independent but the total is modulated by the original 
scintillation variance (equation [4J and so the total residual 
scintillation variance can be modelled by, 

<W = 2a 2 Z0 + ((a 2 scint y x ((a| L ) fc + F 1 )) , (7) 

where cr| L is the scintillation index due to the surface layer, 
F is the Fresnel number used to quantify the 'amount' of 
diffraction and is given by F — D 2 /4Az, and j, k and / are 
solved empirically from the simulation results and are found 
to be j = k = 2/3, I -1.4. 

Using high-resolution generalize d SCIDAR turbu lence 
profile data from San Pedro Martir (jAvila et al.ll2006l ) and 
the model developed from the simulation results we can es- 
timate the expected improvement in intensity variance. The 
SCIDAR profile shown in figure [10] was recorded on 2000 
May 19 and shows a strong turbulent layer at approximately 
10 km throughout the night. Figure [TT] shows the expected 
improvement factor in intensity variance as a function of 
time for the same night. The median improvement ratio is 
11.5 for this example. 

When calculating expected performance for real experi- 
ments it is also necessary to include the exposure time of the 
integration as this will average out the scintillation and re- 
duce the intensity variance. The scintillation index given in 
equation [5] is only valid for very short exposures where there 
is no temporal averaging, i.e. the exposure time has to be 
less than the pupil crossing time of the intensity fluctuations. 
The crossing time, t c , can be calculated as t c = D te i/v w , 
where v w is the velocity of the turbulent layer. If the ex- 
posure time, t, is greater than the crossing time then the 
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Figure 11. Improvement in intensity variance as a function of 
time for the night of 2000 May 19. The median improvement ratio 
for this night is ~11.5. 



scintillation index is modified to jKenvon et alJ l2006). 

" 4/3 rCl(h)h 2 



D, 



^scint (X 



t 



V(h) 



-dh, 



(8) 



where V(h) is the velocity of the turbulent layer at altitude 
h. Using this modification to the scintillation index we can 
calculate an example light curve for a fictional extrasolar 
planet transit for a given turbulence profile. 

Figure [12] shows an example simulated extrasolar planet 
transit. The transit depth is assumed to be 0.05 % and has 
a duration of 2.5 hours. A 2 m telescope and 30 s exposure 
time are also assumed. The optical turbulence profile used 
in the simulation is the same as that shown in figure fTUl as 
measured by SCIDAR at San Pedro Martir. A wind speed 
of 5 ms _1 for the surface layer and 20 ms -1 for all other 
turbulence is assumed. The normalised scintillation noise 
in the visible is reduced from 0.70 x 10 -3 (0.78 mmag) to 
0.21 x 10~ 3 (0.23 mmag), an improvement factor of 3.3. If 
we assume a target magnitude of 11 then we have reduced 
the scintillation to a level which is comparable to the shot 
noise. 

Although the aperture must be placed at the conjugate 
altitude of the turbulence the photometry can be done in the 
focal plane. This means that we do not expect any of the 
other noise sources to increase as a result of implementing 
our conjugate-plane photometry technique. The magnitudes 
of other noise sources, such as shot noise, readout noise or 
flat fielding noise, will depend on other factors. There are 
three possible regimes in which we are interested: scintilla- 
tion dominated, other noise dominated and a mixture of the 
two. In the first and last cases the noise will add in quadra- 
ture and so a reduction in scintillation noise by a factor of n 

will reduce the total noise to, <jt 2 — y ^ o% + ^scint — l) ? 
where gt is the total noise before scintillation reduction. 
Figure [13] shows a 2D plot of the total noise reduction fac- 
tor as a function of the telescope diameter and the target 
magnitude assuming the same parameters as before. The at- 
mospheric model was the median profile from the SCIDAR 
data recorded on 2000 May 19. The optimum telescope size 
is found to be between 1.2 m and 2 m. Less than this and 
the diffraction effects limit the possible scintillation noise re- 
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Figure 12. Simulated light curve of the secondary transit of an 
extrasolar planet with a 0.05% transit depth. The data were cal- 
culated assuming the same atmospheric parameters as measure 
by SCIDAR (figure [TO) and a 2.0 m telescope with 30 s exposure 
times in the v-band and a target magnitude of 11. The top panel 
show the simulated light curves with no scintillation correction 
(black points, top) and with scintillation correction (red points, 
bottom), offset for clarity. The solid lines show the theoretical 
light curve (i.e. with no noise). The data points are randomly se- 
lected from a distribution with a variance equal to the total noise 
at that time, and the error bar indicates the total noise at that 
time. The lower panels show the normalised residuals. 



duction and apertures greater than this become shot noise 
dominated. In the latter scenario the scintillation noise is 
insignificant and so scintillation correction techniques will 
have no effect. 

The median reduction in intensity variance for all avail- 
able SCIDAR data collected over 24 nights in March/ April 
1997 and May 2000 at San Pedro Martir is a factor of 6. 
However, with the limited data available it is difficult to 
say if this representative; it is possible that other times or 
sites will yield even better results if the turbulence is more 
constrained to stratified layers. 

Adaptive optics (AO) can be used to reduce the phase 
aberrations for imaging. Here it is intensity fluctuations 
which are the problem and so AO systems can not directly 
reduce the scintillation. However, AO systems can be used in 
conjunction with this technique to further reduce the scin- 
tillation. As shown previously the surface turbulent layer is 
a major limitation to the conjugate-plane technique. There- 
fore, a ground layer adaptive optics (GLAO) system could 
be used to remove the phase aberrations induced by the tur- 
bulent surface layer and therefore also reduce the residual 
scintillation. On occasions when the atmosphere is domi- 
nated by a number of turbulen t layers a multi-conjugate 
AO system ([Langlois et al.ll2004l ) combined with conjugate 
plane masks could be used to significantly reduce the scin- 
tillation. 



8 J. Osborn et al 




Noise reduction factor Figure 15. Prototype of the conjugate-plane photometer that we 

are due to test on-sky shortly. 



Figure 13. The magnitude of the improvement we can expect to 
observe with conjugate-plane photometry depends on the other 
noise sources. If we assume the same parameters as in figure fT2l 
and an atmosphere given by the median SCIDAR profile then 
we can plot the noise reduction factor as a function of target 
magnitude and telescope diameter. The white line indicates the 
limiting magnitude for a given telescope size to prevent saturating 
a 16-bit analogue to digital converter in a CCD. The optimum 
telescope size is therefore the maximum noise reduction factor 
just above this curve, i.e. between 1.2 m and 2 m. This will vary 
with seeing and camera parameters. 



Conjugate Altitude of 




Figure 14. Conceptual design for one arm of the instrument. 



5 OPTO-MECHANICAL DESIGN 

The design of a conjugate-plane photometer is actually very 
simple. Figure [14] is a diagram of such an instrument in- 
strument. An aperture is placed in the collimated beam at 
the conjugate plane of the turbulent layer. A lens is then 
used to focus the light onto a CCD in the focal plane. As 
the aperture is not in the pupil plane, any off- axis light will 
not illuminate the whole aperture and therefore a separate 
optical arm is required for the target and comparison stars. 
This can be achieved with either a prism near the focal point 
of the telescope, or with pick off mirrors if more stars are 
required. This is completely different to an adaptive optics 
type approach as there are no moving parts once the alti- 
tude has been set. Figure [15] shows the full design of a pro- 
totype instrument, which we shall shortly be commissioning 
to demonstrate the conjugate-plane photometry technique. 



6 CONCLUSIONS 

We have presented a technique, known as conjugate- plane 
photometry, to improve the precision of fast photometry 
from ground based telescopes. The dominant source of noise 
from the Earth's surface is often scintillation due to high 
altitude turbulent layers. By placing an aperture at the con- 
jugate altitude of this layer we can remove the majority of 
the scintillation from this layer. We still detect scintillation 
from other layers, but evidence from turbulence profile mea- 
surements suggests that at premier observing sites the at- 
mosphere typically consists of a single strong high-altitude 
layer and a strong boundary layer. Under such condition our 
technique could remove a large fraction of the scintillation. 
Simulations show that the intensity variance can be reduced 
by an order of magnitude. Theoretical calculations have been 
used to estimate the scintillation noise reduction for a given 
parameter set. For example, with an atmosphere as mea- 
sured by SCIDAR at San Pedro Martir on the 19 th May 
2000, the median reduction in intensity variance is a fac- 
tor 11.5 . Using all available SCIDAR data including times 
when we do not see a dominant high altitude layer we still 
obtain a median improvement of a factor of 6. This is be- 
cause we are reducing the propagation distances from any 
single layer to the conjugate altitude and the scintillation 
index is proportional to propagation distance squared. By 
generating a synthetic light curve for a 2 m telescope in the 
visible using the variance expected from SCIDAR data and 
exposure times of 30 s it was found that we could reduce the 
scintillation noise from 0.78 mmag to 0.21 mmag, compara- 
ble to the shot noise. This reduction in noise will open up 
new science areas from the ground, including the character- 
isation of extrasolar planets through the observations of the 
secondary transit. The conjugate-plane photometer is easy 
to implement as a passive correction technique. However, it 
does require a contemporaneous SCIDAR measurement in 
order to ensure the aperture is at the correct plane. 
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